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Abstract 



^ ' We present an inverse scattering construction of generalised point interactions 

■ (GPI) - point-like objects with non-trivial scattering behaviour. The construction 

is developed for single centre S'-wave GPI models with rational S'-matrices, and 
starts from an integral transform suggested by the scattering data. The theory 
of unitary dilations is then applied to construct a unitary mapping between Pon- 
tryagin spaces which extend the usual position and momentum Hilbert spaces. 
The GPI Hamiltonian is defined as a multiplication operator on the momentum 
Pontryagin space and its free parameters are fixed by a physical locality require- 
ment. We determine the spectral properties and domain of the Hamiltonian in 
general, and construct the resolvent and M0ller wave operators thus verifying that 
the Hamiltonian exhibits the required scattering behaviour. The physical Hilbert 
space is identified. The construction is illustrated by GPI models representing the 
effective range approximation. For negative effective range we recover a known 
class of GPI models, whilst the positive effective range models appear to be new. 
We discuss the interpretation of these models, along with possible extensions to 
our construction. 
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1 Introduction and Main Ideas 



Generalised point interactions (GPI) are solvable models in quantum mechanics rep- 
resenting point objects with non-trivial scattering behaviour. The prototype for such 
models is the class of point interactions (PI), corresponding to Hamiltonians with 6- 
function potentials essentially introduced by Fermi and rigorously defined as self- 
adjoint extensions of —A on the domain of smooth functions compactly supported away 
from the interaction centre @]. This construction leads to 1-parameter families of PI 
Hamiltonians in dimensions 2 and 3 which provide the leading order (scattering length) 
approximation to the scattering behaviour of Schrodinger operators with short range 
potentials in the sector of zero angular momentum (see, for example We refer to 
0] for an extensive bibliography on PI models. 

GPI models are employed to treat more general scattering behaviour, such as higher 
order corrections to S'-wave scattering, non-trivial scattering in non-zero angular mo- 
mentum sectors or point objects in dimensions d > 4. Such models have been studied 
for a long time from the pseudo-potential viewpoint in many body physics - see |[. 
The mathematical study of such generalised point interaction (GPI) models began in the 
Russian literature with the work of Shirokov 0, and more rigorous formulations were 
later developed by Pavlov ||, ^ |10[ and Shondin |[Tl|, [1^ (see also [1^). In contrast to PI 
models, GPI Hamiltonians are not defined on the usual Hilbert space L^(1R'^), but on an 
extended space whose inner product may be indefinite, in which case one must identify 
a physical Hilbert space of states in order to recover the probability interpretation of 
quantum mechanics. 

In this paper, we introduce and develop a new inverse scattering construction for 
single centre GPI models with non-trivial S'-wave scattering. It is currently an open 
problem to extend this construction to higher angular momenta and dimensions d > 4; 
we discuss this further in the Conclusion. Our method is based on the technique of 
unitary dilations due in origin to Sz.-Nagy fl^ and extended by Davis 

The 5- wave inverse scattering problem has been partially studied by Shondin ]TT| , 
as we describe below. First, let us briefly mention the two principal non-inverse GPI 
constructions, which we call the auxiliary space and distributional methods. In the 
auxiliary space method, developed by Pavlov and co-workers 0, |T0[, one starts with 
a given extended Hilbert space and then seeks the class of GPI Hamiltonians which 
'live' on this space. On the other hand, the distributional method of Shondin [0 is 
a direct attempt to define H = —A + |ct;)(co'|, where u is somewhere in the j < 
portion of the scale of Sobolev spaces Hj = (—A + l)~-'/^L^(lR'^). The construction 
leads to a Pontryagin spaceQ H^ = Tio ©C^™, where m is the unique integer such that 
uj e T-C-m-2\H-m-i and the inner product has signature (m, m) on the finite dimensional 
part. The GPI Hamiltonian is then defined on H^ using Krein's formula |[T6[| . 



Returning to the inverse problem, Shondin |TT| considered the inverse scattering 



""^A Pontryagin space is an indefinite (Krein) inner product space with a finite rank of indefiniteness 
see Section 2.1. 
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problem in d = 3 for S-wave scattering data of form 

cot So{k) = k-^r{e), (1.1) 

where r{z) is a rational function with real coefficients. Shondin refers to this class of 
data as the 'i? class': it corresponds exactly with the class of rational S'-wave S-matrices 
satisfying the usual analytic continuation property S{k) = S*{—k) In particular, 

the R class contains truncated low energy expansions of form 

cot 6o{k) = -^ + \kro + kh, + ■ ■ ■ + fc^^+V,, (1.2) 

for low energy parameters L,rQ,ri, . . . ,rn and any n > 0, and thus furnishes approxi- 
mations to the low energy S'-wave behaviour of Schrodinger operators with short range 
potentials to arbitrary order. Shondin's method starts by writing down a candidate 
resolvent on a (positive definite) extension of the usual Hilbert space. (In this respect 
it resembles the auxiliary space method). Various free functions in this resolvent are 
then fixed by requiring that the candidate be the resolvent of a self-adjoint operator. 
However, this method is limited to those r{z) with negative imaginary part in the upper 
half plane, which is a somewhat restrictive sub-class: for example, scattering data of 
form ( p..2| ) is possible only with ri, . . . , r„ = and tq < 0. As we will see later, in the 
context of our method, more general scattering data corresponds to GPI models defined 
on Pontryagin spaces. Thus, in order to apply Shondin's method to such data, one would 
have to guess not only the appropriate extension to the Hilbert space, but also its inner 
product, thereby rendering it much less practical as a construction. 

In contrast, the method proposed here allows one to treat the full R class. It proceeds 
from the simple observation that, for point-like interactions, the scattering data So{k) 
completely specify the S'-wave continuum eigenfunctions as Uk{r) 
So{k)), when scattering normalisation is imposed. This is quite different from the usual 
situation in inverse scattering theory, where the scattering data specifies (in the first 
instance) only the asymptotic form of the Uk{r). As a result, the inverse scattering 
problem for GPI models can be solved without recourse to the usual Gel'fand-Levitan 
machinery. 

The Uk{r) define an integral transform T between Tir = L^((0, oo), dr) and Tik = 
L^((0, oo), dfc), (i.e., the radial position and momentum Hilbert spaces) by 

/ 2 \ ^^"^ r'^ 

{Ti)){k)=^-\ sm{kr + 5o{k))i){r)dr, (1.3) 

We adopt T as a candidate for an eigenfunction transform associated with our desired 
GPI Hamiltonian hgp\. Of course, T cannot be the full eigenfunction transform unless it 
is unitary, in which case we can define hgp\ = T*k^T. However, when T is non-unitary, 
we may 'dilate' it to a unitary operator between enlarged inner product spaces, using 
the theory of unitary dilations as follows. Firstly, we quantify the departure of T from 
unitarity by the operators Mi = U — TT* and M2 = 1 — T*T, and the closures Aii and 
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of their respective ranges. For scattering data in the R class, we will show that the 
rank and signature of Mi and M2 are finite and given in terms of indices reminiscent of 
Levinson's theorem HI? . 



Next, we define indefinite inner products on the Aii, together with an operator 
T : Tir ® M.I ^ Hk (B M.2 which is unitary with respect to the relevant inner products 
and satisfies Pt^^.T'It^^ = T, where Py^f, is the orthoprojector onto Ti^. T is said to be a 
unitary dilation^ of T. We emphasise that the construction of T and the enlarged inner 
product spaces requires no information beyond that encoded in T, and is unique up to 
a unitary equivalence and further dilation. To complete the construction, we define the 
GPI Hamiltonian hgpi by 

where the dagger denotes the Pontryagin space adjoint, and we have used an obvious 
block matrix notation. We will show that A is completely determined by imposing a 
physical locality requirement: that the 'interaction' be localised at the origin. Mathe- 
matically, this is expressed by requiring the Hamiltonian agrees with the free Hamilto- 
nian away from the interaction centre, i.e., hgpi{ip,0)'^ = {—ip",0)'^ if ip E C^(0, 00). 
Subject to this locality condition, we have thus constructed both hgpi and its spectral 
representation. The non-uniqueness in our construction leads to a family of unitarily 
equivalent GPI Hamiltonians with the same spectral and scattering properties. 

Our plan is as follows. In Section 2, we briefly describe some features of analysis 
in indefinite inner product spaces, and also describe the construction of unitary dila- 



tions, essentially following Davis [|T5[. In addition, we sketch our construction in a more 
abstract setting. Next, in Section 3, we explicitly construct the operators Mi and M2 
(which are of finite rank) for a generic subclass of the R class - those whose scattering 
amplitudes exhibit only simple poles on the physical sheet - and compute their rank and 
signature. In Section 4, we construct /igpi as described above. Subject to the locality 
condition, we show that the eigenvalues of /igpi occur at precisely those energies for which 
the scattering amplitude derived from ( |1 . 1| ) exhibits poles on the physical sheet, as is 
the case for ordinary scattering from 'nice' potentials. We construct the corresponding 
eigenfunctions of /igpi, and isolate the physical Hilbert space. We also determine the 
domain and resolvent of /igpi, and explicitly construct its M0ller wave operators in a 
two-space setting, verifying that they exhibit the required scattering theory. 

In Section 5, we illustrate our procedure by constructing GPI models with scattering 
behaviour cot(5o(A;) = —l/{kL) + kM, representing the effective range approximation of 
low energy scattering theory In the case M < 0,Hr and Hk are extended to larger 
Hilbert spaces, and we recover the models of 'type B2' previously constructed by Shondin 



11 . These models also cinSG clS db special case of the auxiliary space construction in M . 
The case M > 0, for which Pontryagin spaces are required, appears to be new. Our 
methods allow the entire class of GPI Hamiltonians to be constructed, along with their 
spectral representations. A particularly interesting subclass of the models constructed 



Sec Section 2 for a note on the nomenclature. 
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corresponds to the case L = oo, with scattering theory cot5o(^) = kM. Such models 
reproduce the leading order behaviour of non-point interactions exhibiting a zero energy 
resonance. We refer to these models as resonance point interactions (RPI). 

We also discuss how these GPI models may be used as models for Schrodinger opera- 
tors with spherically symmetric potentials of compact support. To do this, we employ a 
general methodology for discussing the 'large scale effects of small objects' developed by 
Kay and the author . In particular, we develop fitting formulae (analogous to those 
given in P]) for matching a given potential V{r) to the 'best fit' GPI model. Finally, in 
Section 6, we conclude by discussing various extensions to our method. 

The motivation for the present work arose in a consideration of the scattering of 
charged particles off magnetic flux tubes of small radius [0 , in which it was found that 
the scattering lengths for spin-^ particles generically take the values or cxo in certain 
angular momentum sectors. In consequence, the analogue of PI models representing 
dynamics in the background of an infinitesimally thin wire of flux fails to describe the 
leading order scattering theory in these sectors, and should be replaced by models anal- 
ogous to the RPI models mentioned above. The special nature of this system can be 
attributed to the fact that it is an example of supersymmetric quantum mechanics. 



Elsewhere [19|, we will construct the appropriate class of RPI for this system. 



2 Preliminaries 
2.1 Unitary Dilations 

We begin by describing the unitary dilation theory required in the sequel. Let Hi, ... , Ti^ 
be Hilbert spaces and T G C(Hi, H2). Then T G C(Hi ©7^3, ?^2©'^^4) is called a dilation 
of T if T = P-h^T\t-c-^ where is the orthogonal projector onto H2. In block matrix 
form, T takes form 

Our nomenclature follows that of Halmos ||20|. Elsewhere (e.g., in the work of Davis 
||15|] ) , the term 'dilation' (or 'dilatation') often means that T" is a dilation of T" and 
(T*)" is a dilation of (T*)" for each n = 1, 2, ... (in addition. Hi = 7^2, and Hs = Hi). 
We refer to such operators as power dilations: in the block form ( ^31) , this requires 
PR'Q = for each n = 0, 1, 2, . . .. 

According to a result of Sz.-Nagy [1^, any contraction T from one Hilbert space to 



another (i.e., a bounded operator satisfying ||T|| < 1) has a unitary dilation between 
larger Hilbert spaces. Subsequently, Davis extended this result to arbitrary closed 



densely defined operators at the cost of introducing indefinite inner product spaces. (It 
is clear that if ||T|| > 1, no Hilbert space unitary dilation is possible.) In fact, Davis' 
construction yields a unitary power dilation of the original operator. This has no physical 
relevance in our construction, and so we use a more economical 'cut-down' version of 
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Davis' result, described below. First, we briefly review the salient features of analysis 
in indefinite inner product spaces. Full treatments can be found in the monographs of 
Bognar and Azizov and lokhvidov . 



We employ a particular class of indefinite inner product spaces known as J -spaces. 
Let 7i be a Hilbert space with (positive definite) inner product (■ | ■), equipped with a 
unitary involution, J . We define a non-degenerate indefinite inner product [■, ■] on Ti by 

[x,y] = {x\Jy), (2.2) 

which we call the J -inner product. Ti equipped with the J-inner product is called a 
J-space. Ti, admits decomposition Ti = TC+ © 7i_ = TC+[-\-]TC- into the eigenspaces T-C± 
of J with eigenvalue ±1, where denotes the orthogonal direct sum in the J-inner 
product. If at least one of the TC± is finite dimensional, then TC is a. Pontryagin space 
with respect to [■, ■] . 

The topology of a J-space is determined by the Hilbert space norm; however, operator 
adjoints and the notion of unitarity are defined relative to the J-inner product. Thus if 
Hi {i = 1, 2) are Jj-spaces, and T G £(7ii, 7-^2), the ( Ji, J2)-adjoint of T is defined in 
terms of the Hilbert space adjoint T* by 

= JiTVa. (2.3) 

Equivalently, [T'^x,y]-H^ = [x,Ty]n2 for all x e H2, y e Hi. If [Ux,Uy]'H2 = [x,y]Hi 
for all x,y E V C. Hi, U is said to be (Ji, J2) -isometric; if in addition f/ is a linear 
isomorphism of Hi and H2, and V = Hi, U is said to be (Ji, J2)-unitary. Equivalently, 
UU'^ = l-Hi and WU = 1^2- If T^i = '^a with Ji = J2 = J, terms such as (Ji, J2)- 
isometric are abbreviated to J-isometric etc. 

Returning to the construction of unitary dilations, let T be any bounded operator 
T e C{Hi,H2), and define operators Mi = 1 - TT* and M2 = t - T*T. It is trivial to 
show that the respective closures Aii = RanMj of their ranges are sgn (Mj)-spaces, and 
hence that fCi = Hi® M.i are Jj-spaces, where Jj = © sgn (Mj). We now define a 
dilation T of T by 



^_ , T -sgn(Mi)|Mi|V2 



which has (Ji, J2) -adjoint equal to 



rf^_ jr^j _( T* sgn(M2)|M2|i/2 



(2.5) 



Here, we have used the intertwining relations Tf{T*T) = f{TT*)T and j'*j{j'T*) = 
f {T*T)T* , which hold for any continuous Borel function /. It is now easy to show that 
T^T = l^c^ and TT^ = U/^j) thus verifying that T is a (Ji, J2)-unitary dilation of T. In 
our application. Mi and M2 are finite rank, and so the J-spaces constructed above are 
Pontryagin spaces. 
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We briefly consider the uniqueness of the unitary dilations constructed above. Sup- 
pose Mi are Jj spaces {i = 1,2) and that T : Hi © A/i — > H2 © M is a unitary dilation 
of T with matrix form ( |2.1| ). Then, provided that the Mj are finite rank, one may show 
that 

P-h.^qTU^v = ( J f72 ) ^ ( Ul)' ^^-^^ 

where V = P'I'RanA/i, Q = QRanMi, and Ui and U2 are unitaries (with respect to the 
J- inner products) from JUi and to P and Q respectively. In addition, Ph2®q is an 
orthogonal projection onto 7^2 © Q in 7^2 © A/2. 

Thus T is unique up to further dilation and unitary equivalence of the above form. 
If the Mi are not of finite rank, this statement also holds if the Mj are strictly positive. 
More generally, it is not clear whether Q is necessarily orthocomplemented, and therefore 
whether Pn2®Q exists. 



2.2 Abstract Setting 

In this section, we sketch our construction in a general setting, which makes clear how it 
may be extended. In particular, we show how the domain and action of the Hamiltonian 
is determined. 

Let Hi {i = 1, 2) be Hilbert spaces and let A be a densely defined symmetric operator 
with domain V C Hi. Suppose that A possesses two self-adjoint extensions A± such 
that 

A± = T*AT± (2.7) 

where A is a self-adjoint operator on 7^2 with {A + uj)^^ bounded for some G IR, and 
T± are unitary operators T± : Hi H2- Let a+ and a_ be bounded operators on 7-^2 
which commute with A and define 



r = a+T+ + a_T_. (2.8) 

In our application, a±. are determined by the scattering data. We define Mi and M2 as 
above, for simplicity assuming that they are finite rank (as they are in our application). 
The unitary dilation T derived above is then used to define a self-adjoint operator B on 
the Pontryagin space Hi = 7-^i © A^i by 

where A is a self-adjoint operator on A^2 (with respect to its inner product). Thus 
/ ^ \ / rM(r^-e) + sgnM2|M2|V2A(|M2|i/V + T*U,<l>) \ 

""Ui I -|Mi|V2i(r^_e) + r|^,A(|M2|VV + r*u,$) r 
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where 6 = sgnMi|Mi|^/^$ (considered as an element of 7^2), and B has domain 

D{B) = {(^, $r I - G D{A)}. (2.11) 

To gain a more exphcit description of D{B), we impose the requirement that i? be a 
self-adjoint extension of the non-densely defined operator A©OonD©0 C Hi, i.e., 
B{ip, 0)"^ = {Aip, 0)"^ for all if eV. Later this will carry the physical interpretation of a 
locality condition. It is easy to show that this requirement is satisfied if and only if 
is invariant under A* and 



A = i\M2\-'/'A*\M,\M2\' 



/2y 



(2.12) 



As a consequence of locality, we note that if (y?, $)"^ G D{B) with B{ip,(^)' 



{(f, $)^, then (p G D{A*), and (f = A* (p. For take any ^jj eV. Then 







(2.13) 



J Hi 



We may therefore re-write ( p.ll|) as 



D{B) 



^ G D{A*), 01 G D(i) 



(2.14) 



where 0i = a+7^x+ + '^-'^X- + © and x± = {A±+uj)'^{A* + uj)ip — ip. The advantage 
of this expression is that x± can be shown to be the unique element of ker(yl* + uj) such 
that + x± G D{A±). In our application, x± iiiay be expressed in terms of the value of 
if and its first derivative at the origin. 

To determine the action of B more explicitly, we use the fact that the upper com- 
ponent of the right-hand side of ( 2.101 ) is equal to A*Lp in order to compute = 
sgnMi|Mi|^/2<l>. We obtain 



= -MiA{T^ - 0) + T{A*^ - T*A{r^ - 0)) = TA*^ - A{Tip - 0) 
Using the fact that 0i G D{A), this becomes 

= AQi + w(0i - 0) + r{A* + u)ip - (i + uj){T^ + 01 - 0). 
The last two terms cancel by definition of x± and we conclude that 



B 



(sgn Ml I Ml 1 1/2) -10 



(2.15) 



(2.16) 



(2.17) 



where = y40i + uj{a^T^x+ + ct-^X-)- 
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3 Determination of Mi and M2 



In this section, we determine the operators Mi = 1 — TT* and M2 = 1 — T*T, where 
T is an integral transformation arising from the scattering data in the Shondin R class 



11] given by 



cot (5o(A;) = A;-i4St> H^) =0 for £ > 1, (3.1) 

where p{z) and q{z) are coprime polynomials in lR[z], the ring of polynomials with 
real coefficients. In particular, we will show how the rank and signature of the Mi are 
determined by two 'Levinson indices' defined below. We emphasise that our methods 
are very different to those of Shondin. 

The scattering amplitude corresponding to 5Q{k) is 

Defining the polynomial W{z) by 

W(.\ - / " zq{-z^) p(0) ^0 

^^^^ ~\p{-z^)/z-q{-z^) p(0) = 0, ^"^"^^ 

we note that /o(A;) exhibits poles where = 0. The set Vl of zeros of W{z) in 

the left-hand half-plane Rez < corresponds to poles of ff){k) such that k"^ lies on the 
physical sheet. We refer to the situation where these poles (and hence the corresponding 
zeros of W{z)) are simple as the generic case. In Theorem |4.1| , we will show that the 
discrete spectrum of the GPI Hamiltonian is precisely {E = ~lJ^ | G fi} under the 
requirement of locality.^ 

The qualitative features of the scattering data ( p.l| ) are described by the degrees of 
p and g, two indices defined below, and the asymptotic behaviour of coib^ik) given 

by 

ao = sgn lim cot(5o(fc) and cxoo = sgn lim cot5o(fc)5 (3.4) 

where the limits are allowed to be ±00. The indices are defined by 

J, ^ MO) -Moo) j-,C(0)-CM^ ,3,5j 

TT TT 

where the auxiliary scattering data ({k) is defined as a continuous function on IR^ by 

cotC(fc) = -A:-^|^. (3.6) 

We refer to as the Levinson indices (although Levinson's theorem |]l^ will not hold 
in its usual form). 

■^These eigenvalues can be complex: we will return to this point in section 5.3. 
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We now define the integral transform T = cos6o{k)S + sin5o(^)C, wliicli is suggested 
by tlie naive generalised eigenfunctions Uk{r) = (2/7r)^/^ sin(/cr + 5o(^))- Here, S and C 
are the sine and cosine transforms, defined by 

(S7l))(k) — \ — drtpMsinkr and (Cil>)(k) = \ — drtpM coskr (3.7) 

V TT Jo V TT Jo 

(the integrals are intended as limits in L^-norm). Both are unitary maps from Tir to Tik] 
their inverses have the same form, with r and k exchanged. Thus T is given explicitly 

by 

r- 5 1 c (38) 

{p{k'')^ + k^q{k^)''y/^ {p{k^f + k''q{k'')^y/^ 

Because S and C furnish the spectral representations of —dP/dr^ on L^(1R"'") with Dirich- 
let and Neumann boundary conditions respectively at the origin, we are in the general 
situation of Section 2.2. 

We now restrict to the generic case and explicitly construct the Mj and compute 
their rank and signature. M2 is given by the following proposition, whose proof is given 
later in this section. 



Proposition 3.1 In the generic case, 

= 53a,|U(ed, (3.9) 

where ^^^{r) = e'^^' , and the residue 

— ReSa,22;/o(— i-z). (3.10) 
In addition, RanM2 = spanj^^^ | a; G Q}, and 

rankMs = ^degW^ + 7^ (3.11) 

sigM2 = \(ol-a'0)-ll. (3.12) 

Next, define M.\ to be the space of all L^-vectors of form Q{k^)k{p{k'^)^+k'^q{k'^)'^)~^/^ , 
such that Q{z) G €[z] is a polynomial with complex coefficients. Thus 

Ml = ip{k^y + k^q{k^y)-'^^}^u-i[k^] (3.13) 

where C^f-z] is the r + 1-dimensional complex vector space of polynomials with complex 
coefficients and degree at most r, and 13 — dim Ali is given by 

1 1 

15 = -degW + -{al^- a^) = max{degp, deg g}. (3.14) 

Ml is described by 
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Proposition 3.2 In the generic case, Mi vanishes on A^i"*", and its action on Aii is 
given by MiQ{P)k{p{Pf + Fg(fc2)2)-i/2 = Q(p)k{p{Pf + Pq^Pf)-^!"^, where 

0^e>^Q^eHZ^'-^^^^- (3.15) 

Moreover, Ran Mi = M.i and 

rankMi = ^degW + ^{al - a^) (3.16) 
sigMi = -in + Ij:). (3.17) 

As an example, let us consider the sub-class of the R class considered by Shondin 
11|; namely, the case where r{z) = p{z)/q{z) has negative imaginary part in the upper 
half-plane. In this case, it is easy to show that there can be no solutions to r(— z^) = z 
and hence to W{z) = in the left-hand half-plane, except on the real axis. Moreover, 
one can show that the residues a^^ at these zeros are necessarily positive, so M2 is a 
positive operator as a result of (|3.9|) . Accordingly, T is contractive, and our method 
yields a unitary dilation defined on Hilbert spaces. This explains why Shondin was able 
to construct these GPI models on enlarged Hilbert spaces. 

We now prove the above propositions. 



Proof of Proposition 3.1: M2 may be written in two equivalent forms: 



M2 = S-^ siYi^ 5Q{k)S - C-^ six? 5Q{k)C 

—C^^ sin(5o(fc) cosSo{k)S — S^^ sm6o{k) cos6o{k)C (3.18) 
= C-^ cos^ 6o{k)C-S-^ cos^ 6o{k)S 

—C^^ sin5o(^) cos6o{k)S — S^^ sin5o(fc) cos5o(fc)C (3.19) 

To convert this into an integral kernel we use the following Lemma, which may be proved 
by standard means (cf. Theorem IX. 29 in |]23[). Here, v{x) and w{x) stand for either 
sinx or cosx, and V and W are the corresponding integral transforms from Ti^ to Tik- 

Lemma 3.3 Let g{k) e L2(1R+) n L°°(1R+) and define G = V-^g{k)W. Then G has 
integral kernel 

2 

G{r,r') = - v{kr)w{kr')g{k)dk, (3.20) 



vr JO 

2 



(where the integral is a limit in L -norm). 

In the case degp > degg, sin^5o(^) and sin5o(^) cos5o(^) are fl L°° and so, 
applying Lemma ^]3| to (|3.18| ) and combining terms, M2 has integral kernel 



/-oo ,1 /-oo 7l-/r|'i-2\pifc(r+r') 

M2(r,r') = - / e''^^^'hm6o{k)e'''^^+^'Uk = - ,2, , n 2^ dk. (3.21) 

vr J-00 TT J-00 p[k/) — tkq[k/) 
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Making the substitution z = ik and closing the contour in the left-hand half-plane, the 
integrand has a simple pole at each a; e f2 and (|3.9|) follows. If degg > degp, we argue 
similarly using ( p.l9| ) to obtain the same result as before. 

By linear independence of the and non- vanishing of the a^, it follows that Ran M2 = 
= spanj^j^ I u G Q}, so rankM2 = \ the cardinality of fl. Using residue calculus, 
one may show that 

|n| = ldegM/ + i-rpS*. (3.22) 

' ' 2 ^ 27r i-00 W{ik) ^ ' 

By rewriting the second term as an integral over (0, 00), a small amount of algebra shows 
that the integrand is — 7r~^(5o(fc). Thus is established. 

To compute sigM2, we define the hermitian form m2{ip,ip) : x -^2 — ^ C by 
m2(v?, V^) = {f I M2tp). Labelling the elements of fl as uji, . . . ,w\n\, and writing ip = 
EiQ^o.,, we have 

m2{ip,i)) = Y.'^ii^^^ I I ^^Jcfc = c^E^AEc, (3.23) 

where A and S are hermitian. S has components = {^^. \ ^uj^), and is non-singular 
by linear independence of the C,^. By Sylvester's Law of Inertia 0], the signature of M2 
equals that of A, which has components 

I otherwise. (3.24) 

A has eigenvalues {auj \ u G IR} U {±|a;t^| | u ^ IR}. Labelling the Ui so that ui, . . . ,LJr 
are the real elements of fl, we therefore have sigM2 = sigdiag {a^-^ . . . , a^^). (We have 
used the fact that ajj = a^, and in particular that cj^ G IR implies a,. G IR.) Defining 
C(fc) by (|3.6| ), it is easy to show that cotC(— 1^) = 1 for a; G f2, and that 

= 2 lim N = V (3.25) 

^^-<^ 1 - cotC(^) C (-^) 

Thus sigdiag (ai, . . . , a^) is equal to the number of times that ({k) = 7r/4 (mod vr) as 
k traverses IR"^, counted according to the sign of C'(^) at such points. This is related to 
the Levinson index by ( p.l2| ). ■ 



Proof of Proposition \3. ^ : We compute 



M - cn-i ^^(^ 



(p(F)2 + k'^qik'^Yfi^ {p{k^y + k^q{k^yy/^ 

kq{k^) 1 p{e^ 



{j){k^Y + k^q{k^YYl^^^ {p{k^Y + Fg(P)2)i/2 ' ^^-^^^ 

which vanishes identically on the closure of X' = {p{k'^Y + k'^(l{k'^)Y^'^SC^{Q, 00) as a 
result of elementary properties of the sine and cosine transforms. Furthermore, V is 
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precisely the space A4i defined above, because ip J^V if and only if {piJi^Y + k?' q{k?)Y^'^%lj 
is the sine transform of a distribution supported at the origin and therefore an odd 
polynomial (cf. Theorem V.ll in |^). Hence Mi vanishes on M.^ and Ran Mi C A^i. 

Next, we compute the action of Mi on M.i. By contour integration, 

for polynomials Q{z) such that the operand is in L^. Moreover, it is easy to show that 

from which the action of Mi can be read off as required. 

To compute the rank and signature of Mi, we use the fact that 

rank Mi — rank M2 = sig Mi — sig M2 = dim ker T* — dim ker T, (3.29) 

which follows from the intertwining relations MiT = TM2 and M2T* = T*Mi. It 
therefore remains to determine the dimensions of the relevant kernels. Firstly, note that 
kerT* C Mx and that (from ( ^1271) ) 7/; = Q{k^)k{p{k'^f + fc2g(fc2)2)-i/2 ^ ^erT* if and 



only if ip E Aii and Q^—uj"^) = for each uj E fl. Thus Yiujeni^ + divides Q{z) and 
so 

dimkerT* = min{0 - |(]|,0}. (3.30) 
Now consider ker T. We note that (|3.28|) may be rewritten 



._,2^^. _/2y/2 k p{e)q{-uf) - p{-ul)q{k^) 

Q[ '^J^U (p(p)2 + p^(p)2)l/2 ^2+^2 ' l^-^^J 

and apply the following abstract algebraic result: 

Lemma 3.4 Let Q,R E <C[z] be coprime with maxjdegQ, degi?} = k > 0, and let 
Ai, . . . , Am be distinct elements o/(D. Then the polynomials Pi{z), . . . Pm{z), defined by 

{z - Xi)P^iz) = R{\i)Q{z) - Q{\i)R{z) (3.32) 

span a min{k,m}- dimensional subspace of(Ck-i[z\. 

Proof: Let n = min{k,m}. Then it is enough to show that Pi,...,P„ are linearly 
independent. Assuming that degQ = k, we note that Pi{z) = R{z)Qi{z) — Q{z)Ri{z), 
where Qi{z) = {Q{z) - Q{Xi))/{z - Aj) and Ri{z) = {R{z) - R{Xi))/{z - Aj). Suppose 
the Pi are linearly dependent. Then R{z)S{z) = Q{z)T{z) where S{z) = J^iCaQii^) 
and T{z) = J2i CiiRi{z), for some 7^ (ai, . . . , an)'^ G C". Because Q and R are coprime, 
this implies that S and T vanish identically. But one may easily show that the Qi are 
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linearly independent, by explicitly considering their coefficients. We therefore obtain a 
contradiction. ■ 

In our application, m = \Q\ with Aj = —ujf for each i = 1, ...,m and k = 
max{degp, degg} = U. Thus dimTRanM2 = min{|f2|, 13} and so 

dimkerT = min{|l]| - tJ,0}. (3.33) 

It follows that rank Mi — rank M2 = sig Mi — sig M2 = U — from which (|3.16|) 
and (|3.17|) follow. ■ 



4 The GPI Hamiltonian 

4.1 Locality and Spectral Properties 

The results of the previous two sections allow the construction of a unitary dilation T 
of the integral transform T. Here, we employ T to define a GPI Hamiltonian consistent 
with scattering theory ([31]). We denote H^ = Hr (B Mi and H^ = Hk ® M2 with 
J-inner products specified by J,. = l^^,, © sgn (Mi), and Jk = Iwj. © sgn (M2). In terms 
of our general discussion in Section 2.2, we set A = —dP/dr"^ on domain C^(0, cxd), and 
define = S, T_ = C, setting a+ and a_ to be multiphcation by cosSo{k) and smSo{k) 
respectively. Thus = S*k'^S, the self-adjoint extension of A with Dirichlet boundary 
conditions at the origin, whilst A_ = C*k'^C is the extension with Neumann boundary 
conditions at the origin. The operators A± + 1 both have bounded inverse. 

The S-wave GPI Hamiltonian is defined by 

Vi = ^^( 0' a)^' (4.1) 

where A is a sgn (M2)-self-adjoint operator A"!" = A on A42- To fix A, we require that 
/igpi (■?/', 0)^ = (—■?/'", 0)"^ for all ip ^ ^"^(0, 00) as a locality requirement. For general 
ip G A^2, we have 

= - E I (4.2) 

so is invariant under A* and it follows immediately from Section 2.2 that 
Theorem 4.1 In the generic case, the unique choice of A consistent with locality is 

A = - (sgn(M2)|M2|^/2)"' E ^^^'\UmM2\-'^'. (4.3) 

We proceed to determine the eigenvectors and eigenvalues of A. First note that 
(% I ^2~^^a;,) = ajMii, which follows from the identity = T, ^ul^uj) {^u \ Af2~^^<^»)- It 
is then a matter of computation to see that ipi = (sgn (M2) |M2|^/^) ^^j^ is an eigenvector 
of A with eigenvalue —uf for each i = 1, . . . Because A has rank this exhausts 
the discrete spectrum of /igpi. The following is then immediate. 



13 



Theorem 4.2 In the generic case, and with A. is defined as above, /igpi has the following 
spectral properties: cr{hgp{) = crac(/igpi) U (Xpp(/igpi) where crac(^gpi) = H'*' and <Jpp{hgpi) 
consists of the \ fl\ eigenvalues —ujf, whose corresponding eigenvectors are 

"''' = ^**^-=(r(sg„(Mi|v^)-'e„.)- <") 

The absolutely continuous subspace is the Hilbert space T^Tik- 

This bears out our earlier statement that the poles of the scattering amplitude on 
the physical sheet correspond to the discrete energy spectrum, if locality is imposed. 

The physical Hilbert space is required to be a positive definite invariant subspace of 
IIj. relative to /igpi.0 In 11^, we have the [■, -jnfe-orthogonal decomposition 11^ = 'Hk[+]M.2i 
where is spanned by the eigenvectors (fi of A. We compute 

= (U I M,%^) = I ^ ^ ^ (4.5) 

Hence H^ is decomposable as H^ = 7Yfc[+]-^+[+]-^-[+]-f^ where is spanned by the (pi 
with [ipi, ipi]M2 > > 0)) ^- is spanned by those with [ipi, ^i\M2 < ('^o;; < 0), and 
H is the hyperbolic invariant subspace spanned by those with uji IR. Moreover, this 
is a decomposition into invariant subspaces, because D{k'^) is dense in Tik- The physical 
Hilbert space Hphys is therefore defined by 

^phys = fHnk[+]E+). (4.6) 



We briefly discuss the uniqueness of the GPI Hamiltonian constructed in this way. As 
noted in Section 2.1, T is unique up to further unitary dilation and unitary equivalence 
because the Mj are of finite rank. Further dilation merely corresponds to the (trivial) 
freedom to form the direct sum of /igpi with the Hamiltonian of an arbitrary independent 
system. On the other hand, replacing T by (1 © U2)T{t © Ui) where Ui is a sgnMj- 
unitary operator on Aii for i = 1,2, it is easy to show that the local GPI Hamiltonian 
/igpj obtained is given by 




We have therefore constructed a family of unitarily equivalent GPI Hamiltonians on H^ 
corresponding to the same scattering data. It is clearly sufficient to study hgpi alone in 
order to determine the domain and scattering properties of h'^^^. 

*An invariant subspace £ of a J-space IC relative to a linear operator A on is a subspace of /C 
such that D{A) D C = C and RanAj^ C C, where the closure is taken in the norm topology of IC. 
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4.2 Domain and Resolvent 



We now determine the domain and explicit action of the operator hgpi under the locahty 
assumption. Our result is the following: 

Theorem 4.3 Let Gq = {2/7iy^^k'^^-^{p{k'^f + k^q{k^f)-'^/'^ . Then in the generic case, 
D{hgpi) - I (| ^ ) k> V'' e ACioc(0, oo), ip, ip" e L^; $ e Ml, 

e-A[(^]eo e L'(A;')nMi|, (4.8) 

where © = sgnMi|Mi|^/^$ and 



AM 



P(p{0) degp>degq 
Pip{0)-Qip'{0) degp^degq (4.9) 
-Qip'iO) degp<degg, 



and P and Q are the leading coefficients of p{z) and q{z) respectively. (In the case 
Ml = 0, D{hgpi) = {</? I </?, </?' e ACiociO, oo), (f, (p" e L^; X[(p] = 0}.) Moreover, 

"SP' ( S ) - ( 1" ) . ("0) 

where $ is given in terms of Q — sgnMi|Mi|^/^$ by 

e = Me-AMeo) + y (^PFT^MW^^ ' ^ ^ 

Proof: The result is a direct application of the discussion in Section 2.2. The key point 
is that, for each ip G D{—d^/dr^\^^^Q^^), the vectors x+ ^^'^ X- given by 

X+^-me-' and X- = ¥''(0)e-^ (4.12) 

which follows because x+ (X-) is the unique element of ker(— c/^/dr^l^oo^g oo) + 1) such 
that (p + x+ i'P + X-) is in the domain of the Laplacian with Dirichlet (Neumann) 
boundary conditions at the origin. ■ 

The resolvent of /igpi may be written in the form of Krein's formula as 

(V. - - ( ' ) + ,(,),g)v.,(,) ^W^W* ("3) 

Here, Ro{z) — S'^ik"^ — z)~^S is the free resolvent and the defect element F{z) e 11^ is 
given by 

^^^^ " [ (sgnMi|Mi|V2)-i^(^) ) ' (^-1^) 
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where "^(z) E Aii is 



2^/2 k{p{e)q{z) - p{z)q{e)) 



(P-;z)(p(A;2)2 + A;2g(A;2)2)i/2' 
and the operator Ri{z) is defined on A^i by 



(4.15) 



(4.16) 



where Q(2;) is defined in terms of $ by 
e = sgnMi|Mi|^/2$ = 



2. 1/2 



TT 



kQ{e 



(4.17) 



The above expression for R{z) may be verified directly using Theorem [4.3| , and the fact 
that 

[(sgnMi|Mi|V2)-ivI;(^),$]^^ = ^^^^g) 

which is required when one takes inner products with F{z). Using this result, it follows 
that ( [4.13|) holds for elements of form (0, $)^ with Q{z) = 0; direct computation estab- 
lishes it for Q{z) = 1 and also for vectors of form {ip,0)'^ with ip G 7^r-Q Thus ( [4.13|) 
holds on the whole of 11^. It remains to establish equation (^4.18|) . Multiplying through 
by q{z), the LHS of ( [4.18| ) is equal to 



(g(z)^(z) I A'/f ^0) = (r*g(^)^(z) | M-'r*Q) + (g(z)^(z) | 6) 



(4.19) 



Using the identity \ M2 ^^uji) = Ci^Hij and the results of Section 3, the first term is 

(4.20) 



{T*q{-z)^>{-z) I M,-^r*e) = E 



^eQg(-w^)(w^ + 2: 
The required result then follows from the calculation 



{q{-z)^{z) I 6) 



1 /"^ ikQ{P){p{z) - ikq{z)) 
{k^ - z){p{k^) - ikq{k^)) 



TC J- 



\ P{z)-ujq{z) 2\ 



g(— C<j2)(ct;2 _)_ 2;) 



(4.21) 



4.3 Scattering Theory 

In this section, we construct M0ller wave operators for hgpi relative to the free Hamil- 
tonian ho = S^^k^S on Ti^ in order to check that /igpi actually exhibits the required 
scattering behaviour. Because scattering is a function of the continuous spectrum only, 

^Here, it is useful to employ the decomposition Tir — Ran (— d^/fir^ — 2;)|cj=(o,oo) ® Ce"^"^^^^^''. 
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our results in this section are actually independent of the precise form of A, and therefore 
of the locality requirement. 

We work in the S'-wave, and employ a two space setting: let B be self-adjoint on Hi, 
A be self-adjoint on H2 and J" be a bounded operator from Hi to H2- Then the M0ller 
operators fl'^^A, B; J') are defined by 

Q^{A,B-J)= lim e'^'je-'^^'P^B), (4.22) 

t — ^=poo 

and are said to be complete if the closure of Ranfl'^^A, B; J') is equal to RanPac(^)- 

In the following, J'r and J'k are the natural embeddings of Hr and Hk into 11^ and 
Uk respectively. 

Theorem 4.4 Let J : Hr — > H,. he given by J = 'f'^ J^T ■ Then ^'^{h^^\, ho; J') exist, 
are complete, and given by 

^^{h^pi, ho; J) = f^Jke^''''-''^S, (4.23) 
where So{k) is given by ( ^7^) . 

Proof: Writing Ut for multiplication by e^*^^* on Hk, we have 

= f^JkU-tTS-'UtS. (4.24) 

Now, for any u{k) G C^(0, 00), 

\\U_tTS-^Utuik) - e±'^"(%(A;)f = || smSoik)C{C-^ ± tS^^)Utuik)f 



2 r'^ 
< - dr 



TC Jo 



dke'^^^'-^''\{k) 



(4.25) 



which vanishes as t — =fcxd by (non)-stationary phase arguments (see the Corollary to 
Theorem XI.14 in [^). Thus U^t'^S^^Ut Q±^Soik) g^rongly as t — > =Foo. The existence 
and form of the M0ller operators are then immediate. One easily checks that they are 
unitary maps from Hr to Pac(^gpi) = JkHk, to establish completeness. ■ 

We conclude that our construction does indeed yield the required scattering theory, 
and also that - as a by-product of the construction - complete M0ller operators may 
easily and explicitly be determined. 



5 Examples 

As an application, we construct the class of GPI models with scattering data 

coi5o{k) = -^ + kM, (5.1) 
kL 
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where L is the scattering length, and M is twice the effective range. These models 
therefore represent the effective range approximation to the behaviour of a non-point in- 



teraction in the S'-wave. This class of models has been partially studied by Shondin |TT 
who considered the case M < ('models of type -B2') and also appears as a special case 
of the models considered by Pavlov in (We also note that van Diejen and Tip 
have constructed models of type cot6o{k) = {ak + bk^ + ck^)~^ using the distributional 
method.) The case M > does not appear to have been treated before. Our construc- 
tion provides a unified construction for all models in the above class, and also provides 
the spectral representation such models as a by-product of the construction (although 
we will not state this explicitly). 

The above class of GPI models contains two interesting sub-families: the ordinary 
point interactions (M = 0) and also the resonance point interactions arising formally by 
setting L = 00, i.e., cot(5o(/c) = kM with M G IR U {00}. Such models are required in 
situations where the scattering length is generically forced to be infinite, for example in 
certain systems of supersymmetric quantum mechanics. 

We begin by briefly treating the point interactions, both for completeness and also 
to demonstrate how this class arises in our formalism. We then turn to the general case, 
obtaining RPI models in the limit L ^ —00. 

5.1 Point Interactions 

The required integral transform is 

r = (1 + {kLfy^/^s - kL{i + {kLfy^i^c. (5.2) 

In the cases L = 0, cxo, T reduces to S and C respectively, and the Hamiltonian is given 
immediately by T*k'^T. We exclude these cases from the rest of our discussion. 

We therefore apply the construction of Section 3, with p{z) = —L^^ and q{z) = 1. 
We find that L5 = 0, so Mi = (i.e., TT* = 1). Straightforward apphcation of 
Proposition yields 

M - I \Xl){Xl\ L>Q 

M2-j Q ^^Q^ (5.3) 

where X-l('") = (2/ LY^'^e~^^^ is normalised to unity. Hence if L < 0, T is unitary and 
the Hamiltonian is hi = T*k'^T, with purely absolutely continuous spectrum IR"*". In 
the case L > 0, the momentum Hilbert space is extended to Hk © C, representing a 
single bound state, and the unitary dilation T : Hr ^ 'Hk © C takes form 



^-itelj^ t.^(r. |x.>). (5.4) 

(Tifc ©C has the obvious inner product.) The Hamiltonian is 

hL = f-'(^' \ \f = T*eT + \\xL){XLl (5.5) 
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and the locality requirement fixes A = — which is, of course, the usual value. Finally, 
the domain of /i^ is given by Theorem ^]3| as the space of 99 with G ACiodO, 00), 
ifi" G and satisfying the well known boundary condition 

¥?(0) + V(0) = 0. (5.6) 

To summarise, all the well known properties of point interactions may be derived within 
our formalism. 



5.2 Effective Range Approximation 

In this section, we maintain M 7^ 0, L 7^ 0, setting p{z) = —L~^ + zM and q{z) = 1. 
We will not explicitly construct the dilation (although this follows immediately from our 
discussion), but will use the results of Section 4 to read off the domain and action of the 
GPI Hamiltonian h^^M- 

Using the results of Section 3, we find 

= 1; N = | ^ + !lf"^ + r^^ f """^ (5-7) 
' ' [ ^(1 + sgnM) L = 00. ^ ^ 

Writing W{z) = —M{z — uji){z — UJ2), ft is the subset of {ci;i,u;2} lying in the left-hand 
half-plane, and we have ui + uj2 = —M~^, ujiuj2 = {ML)~^. The residues a^^ are 

"^2 = T77-. — T- (5- 



In addition, the space A4i = Ran Mi is equal to €\r]), where 



^(^)=^(F + (FM^-L-T)^/^' ^'-'^ 



and the normalisation constant is 



Af- (2|M|/7r)V2 ML>0 

~ 1 (2|M|/7r)V2(i_4ML-i)V4 ML < 0. 



Using Proposition 3.2, we obtain 



+1 M<0,L<0 
Mi = \\r]){r]\; X={ -sgn M(l - 4ML-i)-V2 ML < (5.11) 

-1 M>0,L>0. 

Accordingly, the extended position inner product space is 11^ = Tir © C with J-inner 
product specified by J = 1 © (— sgnM). The scalar component is the coefficient of \ri) 
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in Ail. For all generic cases (i.e., all cases other than L = AM > 0) Theorem f4.3| entails 
that the domain of hi^M is 

^(^l,m) = |(^ I ^ \^,^' eAQociO,^), ^,^"eL^; $ = -|M|i/V(0)|, (5.12) 
and that the action is 



J' 



^l,m[ _|m|V2^(o) ) - ( -sgnM|Mri/2(^'(0) + L-V(0)) ^^'^^^ 



Moreover, one may show that these equations also hold in the non-generic case L = 
AM > 0. 

It is worth noting how this domain and action correspond to the scattering data ( p.l[ ). 
Solving the equation /il,a/(v'5 = k'^if, ^)'^ for the generalised eigenfunctions of h^^M, 
we find (p{r) oc sin(A;r + d{k)) for some d{k), and also obtain the relation 

- sgnM\M\-^/^{^'{0) + L-V(O)) = ~k^\M\'/^^{0), (5.14) 

which entails that k cot d{k) = ip'{0)/Lp{0) = —L^^ + k'^M. Thus d{k) is precisely the 
scattering data d^ik). 

The RPI models, which have scattering data cot(5o(A;) = kM are obtained in the 
same way. The space Mi is spanned by ipM{k) = (2|M|/7r)^/^(l + (A;M)^)~^/^, and the 
operator Mi is found to be Mi = —{sgnM)\tpM){4'M\- Thus the inner product space is 
Hr = Ti.r ©C with J = 1 © (— sgnM). They have the domain ( p.l2| ) and action 

( -|M|VV(0) ) = ( -sgnM|M|-VV(0) ) ' ^^'^^^ 

Let us consider the physical Hilbert space for these models. From Section 4, this is 
constructed by projecting out the hyperbolic invariant subspace, and also those eigen- 
functions with negative norm squared (if present). The bound states of h^ M are clearly 
vectors of form (^^, |M|^/^)'^ with norm squared equal to — (2Reco')~-'^ — M, where cu is a 
root of uj"^ + M~'^uj + {ML)~^ = 0. There are four cases to consider: 

Case (i): M < 0. 11^ is positive definite so no projection is required. 

Case (a): M > 0, L < 0. There is a unique bound state with 

l^Ml-^^ (5.16) 

-2M ^ ' 



and negative norm squared. Projecting this state out, we obtain 



^phys= . (5.17) 
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Case (Hi): M > 0, < L < AM. There are two bound states with complex conjugate 
eigenvalues. Accordingly, their eigenfunctions span a hyperbolic invariant subspace. 
Projecting this subspace out, we find 

Hpv = { ( ^ I ) I ^ e H, s.t. - &) | ^) = o) . (5.18) 



where uj is given by equation ( |5.16| ). 



Case (iv): M > 0, L > AM. There are two bound states with real eigenvalues. However, 
only the state specified by ( [5.16| ) has negative norm. Projecting this out, we arrive at 



the same expression for Tiphys as in case (ii). 

RPI models are covered by Case (i) for M < 0, and have 7-^phys given by ( p. 171) for 
M > 0, with u = -1/M. 

The GPI Hamiltonian acts on Tiphys by restriction. For example, in case (ii) above, 
we have 

D{hL,M\Hp^y.) = I ( I y?) ) ' ^' ^ ^<^ioc(0, oo), ip, if" e L^; 

M<^(0) = -(ej^)| (5.19) 

on which /iL,Af Iwphys ^"^^^ before. The restricted operator has the same continuum 
spectrum as Hl^m, but has no bound states in this case. Moreover, the property of 
locality is partially lost: it is clear that vectors of form (99, 0)"^ with ip G C^(0, 00) are 
in Tiphys only if ± However, for elements of this form in Tiphys, it remains the case 
that /iL,A/|wpi^y^(v5, 0)"^ = {—{p",0)'^. Thus the properties of locality and 'positivity' are 
not entirely compatible. 



5.3 Physical Interpretation 

In this section, we discuss how the effective range models constructed above may be used 
to model Schrodinger operators H = —A + V, where V is smooth, spherically symmetric 
and compactly supported within radius a of the origin. Our methodology extends that 
described in [^], in which the scattering length approximation is discussed. 

Given a smooth spherically symmetric potential V{r) supported within radius a of 
the origin, we may find the 'best fit' GPI model /il.a/ as follows. Let uq be the S'-wave 
zero energy eigenfunction, i.e., the solution to —Uq + Vuq = with regular boundary 
conditions at the origin. Then the arguments of Section 11.2 of give the low energy 
parameters L and M as 

L = a 



(5.20) 
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and 

Thus the scattering behaviour is cot5o(^) = —{kL)^^ + kM + 0{k^) and the best fit GPI 
model in our class is hL,M- We refer to equations ( ^.20| ) and ( |5.21| ) as fitting formulae; 
equation ( ^.20|) is the fitting formula employed in 0. The range of energies for which 



the approximation is valid can be determined by a 'believability' analysis analogous to 
that described in Q. We will not do this here. 

Note that M obeys the bound 

-oo<M <a\l-^ + l (^] I . (5.22) 



Moreover, this bound is best possible: for any L G IR U {oo} and any M in the above 
range, one can clearly find a smooth function UQ{r) satisfying regular boundary condi- 
tions at the origin, uq oc (1 — r/L) for r > a and such that (|5.21| ) holds. Then the 



potential defined by V{r) = UQ{r)/uQ{r) has S'-wave scattering behaviour approximated 
to second order by hL,M- The contribution to the total scattering cross section from the 
effective range term generally outweighs that from higher angular momenta, so the S'- 
wave GPI model provides a second order approximation to the full scattering behaviour. 

Finally, we discuss the interpretation of the discrete spectrum of h^^M- We have 
constructed h^^M so that its scattering behaviour matches that of a given Schrodinger 
operator at low energies, E. For larger \E\, the approximation breaks down - in the 
language of P we say that it is no longer 'believable'. Thus, deeply bound states are 
unlikely to be believable. In particular, for < L < 4M, hi^M exhibits a complex 
conjugate pair of eigenvalues, which can never be believable.^ Such phenomena are 
artifacts of the idealisation process, due to the truncation of the low energy expansion. 
The issue of believability is discussed in [Bl; similar comments are made in |13 . 



6 Conclusion 

We begin by discussing various generalisations of our method. There are many situations 
in which the analysis of Section 2.2 may be applied. In two dimensions, for example, 
one can consider radial GPI Hamiltonians which agree with 



h = —-r- + - (6.1) 
r dr dr 



away from the origin, as models for an infinitesimal 'dot' of magnetic flux z/, with |z/| < 1. 
In this case, one must employ Hankel transforms rather than sine and cosine transforms. 



^These are not eigenvalues of hL,M restricted to the physical Hilbert space. However, they per- 
sist as poles in the scattering amplitude and our remarks still apply: Hl^aj does not give a reliable 
approximation to the scattering theory at those scales. 
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In [|I9| we will implement this programme to construct a models of RPI type for the 
Dirac equation in the presence of an infinitesimal tube of flux. These models provide 
the leading order approximation to the scattering data. 

Our method could also be applied to S-wave GPI models with a Coulombic tail. In 
this case, the appropriate integral transforms would be based upon Whittaker functions 
and the scattering data would be specified in terms of Coulomb-modified partial wave 
shifts. In this case, the dimension of A^2 would be countably infinite, due to the count- 
able discrete spectrum of such models. However, one would expect A^i to remain finite 
dimensional for simple models. 

Secondly, it is of interest to generalise the unitary dilation method to sectors of higher 
angular momentum with £ > 1 (and the corresponding analogues for magnetic flux dots 
- i.e., |z^| > 1 - and Coulombic GPI for ^ > 1). This is more problematic, because the 
radial Hamiltonian —d'^/dr'^ + i{i+ l)/r^ is essentially self-adjoint on C^(0, oo) and so 
the method of Section 2.2 does not apply. Here, it might be possible to obtain a suitable 
integral transform by analysing the distributional construction. We hope to return to 
this elsewhere. 

Finally, we consider applications to the definition of arrays of point scatterers. Here, 
the most likely use of our methods is to generate the 'monomer' by inverse scattering. 
By passing to the resolvent written in the form of Krein's formula, one can isolate the 



appropriate 'defect element' and proceed to form the array by methods discussed in |]T3 
which generalise the procedure for arrays of PI developed in p7|. 



To summarise, we have introduced an inverse scattering construction for GPI mod- 
els using the theory of unitary dilations, and developed the method in detail for the 
class of single centre S'-wave GPI models with rational ^-matrices. A physical locality 
requirement completes the specification of the Hamiltonian, whose scattering, spectral 
and domain properties are explicitly determined from our results. 
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